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ABSTRACT 

The antiplane shear problem for two bonded dissimilar half 
planes containing a semi -inf ini te crack or two arbitrarily located 
collinear cracks is considered. For the semi -i nfi ni te crack the 
problem is solved for a concentrated wedge load and the stress 
intensity factor and the angular distribution of stresses are 
calculated. For finite cracks the problem is reduced to a pair 
of integral equations. Numerical results are obtained for cracks 
fully imbedded in a homogeneous medium, one crack tip touching 
the interface, and a crack crossing the interface for various crack 
angles. 
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1. INTRODUCTION 

In analytical studies relating to the fracture of structural 
solids, solutions of anit-plane sfiear probTetns serve two distinct 
purposes. First, they may be: used to shed some light on the quali- 
tative behavior of the solutions for the corresponding somewhat 
more difficult in-plane deformation problems. Secondly, they may 
have practical applications in their own right in situations such 
as torsion or three-dimensional problems in which the so-called 
third mode is separable. Particularly in the fracture of nonhomo- 
geneous materials initiating from the flaws in the close neighbor- 
hood of phase boundaries, such as in the growth of microflaws in 
polycrystals and cracks in composites, the third mode of fracture 
may be quite important. With this in mind, it seemed to be worth- 
while to report some new results regarding the anti-plane shear 
cracks in nonhomogeneous materials. The main problem is that of 
a slanted crack in the neighborhood of a plane bimaterial interface, 
including the geometries of the crack touching or crossing the 
interface. The cases of semi -i nf i ni te crack and the half plane 
consisting of two bonded wedges are also considered as they relate 
to the main problem. 

The general antiplane shear problem is wel 1 -understood and the 
basic techniques dealing with its solution have been thoroughly 
documented in literature. Therefore, to include in this paper an 
exhaustive review of the subject and the details of the analysis 
did not seem to be necessary. The problem under consideration is 
the generalization of that considered in [1] where only the case 
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of the cracks perpendicular to the interface was studied. The 
general method has also been described in [2]. References on the 
crack and dislocation problems may be; found in [1] and [2j. 


2. SEMI-INFINITE CRACK AND THE HALF PLANE 

The anti-plane shear problems for bonded dissimilar half 
planes with a semi -inf inite crack and for a half plane which con- 
sists of two dissimilar wedges are described in Figures 1 and 2. 

In the case of semi -i nf i n i te crack under a concentrated anti-plane 
shear loading given by 




0 ) 


using the standard Mellin transforms the solution may be found as 


w^(r,e) - 5^ 


rc + i“ 


[A.(s)sin(s0)+B.(s)cos(s0)]r"^ds, 

c - i 00 -I J 


■ 2Vi 


c+i“ 


C-l“> 


vi.s[A.(s)cos(s0)-B.(s)sin(s0)]r 

J J J 


-s-1 


ds , 




f c + i 00 


c- V 


vi.s£A.(s)sin(s0)+B.(s)cos(s0)Jr"^'^ds, 

J J J 


(j = 1,2,3) . (2a-c) 


where w. is the z-component of the displacement, x. and x. are 
J J02 Jy.2 

the stresses in the jth wedge shown in Figure 1, and 


' "s'') ' , 

-y|A.| (s ) s i n. (2sTr) 


A^Cs) = 


s i n ( s TT ) K ( s ) 
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B2(s) 


(s)[1-cos(2sit)] 


Bo(s) = 


s i n ( s ir ) K ( s ) ’ 

A2(s )sin[s(0Q+Tr)] + B2(s )cos[s(0Q+Tr)]-A3(s)sin[s(eQ-Tr)] 


cos [s ( 0 - -tt) I 


B,(s) = 


A2(s)sin(s0^)+B2(s)cos(s0Q)-A^ (s )sin(s0Q) 


cos(s0^) 
0 


(3a-e) 


K(s) = (y^ 


(4) 


Of particular practical interest is the shear cleavage stress in 
the half plane 2 which, for r<<r^, may be expressed in the following 
asymptotic form: 


T, (r.e) = ^ 

0Z 0 


G,-|(e) ^ 


0(r/rQj^2‘^) , B=l-s^ 

0o<0<6o + 'n' 
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r<<r^ 

0 


> 


(5) 


-2y2COs[s^ (iT-9)] 

" (y ^ -y 2 ) ( TT- 20 Q ) s i n Ls ^ { tt- 2 ) J - (p l+y 2 ) irs i n ( s ^ n) ’ 

(0o-®-®o''’^^ * (6) 

where s is the root of K(s) = 0 with 0<Re(s-|)<l and it can be shown 
that there is only one such root in this strip which is always real, 
and Re(s 2 )>l. Table 1 shows the power of stress singularity for 
various material combinations and for various values of 0^. Expres- 
sions similar to (5) can be developed for other stress components 
in the two materials. Figures 3-5 show some sample results for the 
function G ^(0) in the entire range O<0<2it. Analytically, it can 
be shown that G ^(0) becomes maximum for 0 = tt indicating that the 
plane of the crack is the weak shear cleavage plane in the bonded 
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medium. One may also note that (arbitrarily) defining a stress 
intensity factor by 

^ ^ 20Z 

it is seen that in this problem 

S " q ''o’^ '^-1^''^ • 

Also, referring to Figure 1 and defining the stress intensity factor 
in terms of the crack opening displacement as (see [1]), 

•^3 = ^*‘9r (5) 


we find 

I* 


sinLs^ (7T-e^)JcosLs^ (7r-0^)] + sin(s^ 6 q)cos(s^ 0q) 


( 10 ) 


Figure 6 shows an example for the angular distribution of the 

asymptotic values of which, for the nonhomogeneous medium, is 

seen to be discontinuous at the interfaces 0=0 and 0 = rr + 0 . 

0 0 

For the half plane shown in Figure 2 the basic equations (2) 

are still valid where for the external load 

^iez<''’-®o> ‘ qi«(r-r„). (11) 

the functions A^.(s) and B^.(s), (i = l,2) are given by 

q,r^ sin(s0.) p 
A,(s) = ^ A,(S) . 




‘1 


q,r^ cos(s0.) 

®l<^) = - u^sT T f T r" = ^^(s) 


T^sXjTsJ 


0p = TT-0^ 


( 12 ) 


Ki(s) = 2 ^{ (y^+P 2 )sin(sTT) + (y^-y^ ) s i n [s (tt-20^ ) ] } 


(13) 
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Table 1. Power of stress si ngul ari ty /3 , for a crack 
terminating at the interface. ” 0.0072 for 

Boron-Epoxy, 0.043 for Aluminum-Epoxy, 1.0 for homo- 
geneous medium, 23.08 for Epoxy- A1 umi num , and f38 . 46 
for Epoxy-Boron. 



Table 2. Power ofsingularity 3' at vertex 
for bonded wedges of total angle tt. 
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Examining the asymptotic behavior of the solution for it 

can be shown that for e^<n/2 the power 3' of the dominant term is 
negative if M-j<P 2 (i-®*> stress state at r = 0 is bounded), 

and positive if stress state at r = 0 is singular). 

Near the apex the shear cleavage stress is (Figure 2) 


(r,e) 

lez 



“i(®) s 1 

+ 0(r/r ' ) 

(r/r “ 
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H^O) 


(y^-y2)sin[s.|(92-e)]-(y.|+y2)sin[s^(e2-0)] 



(14) 


-2y.sin[s, (e.,-0)] 

«2'®) = ' - v ri^ T ^ — 


0o = tt-9 

C 0 


(15) 


L(s^ )=(y^+y2)7TCOS (s^ 7r) + (y2-y^ ) (7r-2e^)cos[s^ (tt-20q)] , (16) 

It can also be shown that the displacement derivatives and the 
cleavage shear along the bond line 0=0 are related by 

^1 Jr ^ {(y^ -y2)cos[s^ (tt-20^)] 

+ (y^+y2)cos s^tt} W 2 (0 ,7t- 6q ) . (17) 


ii"ir->o'"^ 'Ti (»",0) - lim^^Qy^ r^ 
y z 


^1 " (y 2 ~y^ )co 


-2y^y2Sin[s-| (tt-6^)] 


i(tt-20^, 


,y^+y2)cosTs^TT 


Table 2 shows some sample results for 3' where the modulus ratio 


used as examples correspond to Aluminum-Epoxy and Boron-Epoxy 


material pairs. 
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3. FINITE CRACKS 


Referring now to Figure 7 we consider the bonded half planes 
containing two collinear cracks (a-]<r<b^ , 0= Q-j ) and (a 2 <r<b 2 , 6=tt+9i ) . 
It will be assumed that the crack surface tractions 


~ ^29 ( 20 ) 
are the only external loads acting on the medium. The results for 
all other types of loading may be obtained through the standard 
superposition technique. Defining the unknown functions 

fi(r) = ^[w.(r,e.+0) - w.(r,e.-0)] . 02 = tt + 9^ , 


i=l,2 (21) 

noting that (see Figure 7) 

f^. (r) = 0 , 0<_r<a^ , b^<r<“ , i = l,2 (22) 

and using the standard dislocation solutions as the Green's func- 
tions [3], the problem may be reduced to the following system of 
integral equations: 


^l^P^ ^ j ^ ^l^^^^s^ ■ k(p,s,0^ )}ds 


- ^ 1 / ^ 2 ^"^ ^ ’ (a^<P<b^) 


^2^*^^ ^ 2 I I ^ ^2^^^ ^s^ k(s,p,0.| )}ds 


- (1-X) 


U2 


^ ^l(s) ^ ds, (a 2 <P<b 2 ) (23a, b) 
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where 




( 24 ) 


,,, . „ , , X 

k(p,s,6i) . ,o-S\5 ? 9 

(^)^sin20i + cos^e^ 


(25) 


From (21) and (22) it follows that (23) must be solved under the 
following s i ngl e- val uedness conditions: 

-b. 
a. 

The solution of the system of singular integral equations is 
of the form 


/•t>i 

J f^ (r)dr = 0 , i = l.,2 . 


(26) 


1 


G.(s) 

f.(s) = ^ — 

^ (b .-s)°‘j (s-a 


, j = 1.2 


(27) 


where G.(s) is a bounded function in the closed interval a-<s£b.. 

J J J 

The coefficients a. and $. may be determined from the application 

J J 

of the function-theoretic method to the singular integral equa- 
tions [4]. Thus it can be shown that 

6i = 32 = ]■ > for a^>0 , a2>0 , 

3^ = 3 , 62 = |- . for a^ = 0 , a2>0 , 

3^ = 32 = 3' , for a^ = 0 = a2 » (28a-d) 


where the real numbers 3 (0<3<!) and 3' (0<3'<0.5) are the powers 
of stress singularity found in the previous section for a semi- 
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infinite crack and for a nonhomogeneous half plane, respectively . 
(see tables 1 and 2). 

The integral equations (23) are solved by introducing the 
following normalized quantities 




2s - (b. + aj 


2p - (b . + a . ) 
3 y . _ ^ Vl 3 ^ 

b,. - a,. * “ 


3 3 


■ ^3 


fj{s) = 4.j(t) = 


9^ 


b,-a. “j+®j “j 

(^V^) (1-tj) (1 + tp 


qj(p) = Qj(x) « »tj.)<i ) , (j = 1,2) 


(29) 


and by using Gauss-Chebyshev or Gauss-Jacobi integration formulas 
(see, for example [5]). For crack geometries a^>0 and a^^O or 
a-]^0 , a 2>0 the conditions (26) are valid and the solution is 

rather straightforward. For a-j = 0 = a 2 (26) reduces to the fol- 
lowing condition 


f, (s)ds + f 2 (s)ds = 0 

J n J n ^ 


(30) 


Stating that the crack surface displacements at r = 0 must be 
continuous. The additional condition which is necessary to insure 
the uniqueness of the solution of (23) comes from the function 
theoretic analysis (see, for example [5]) and is given by 


P]G^(0) = -G2(0)/bpF2j(y^-y2)cos[2(l-B^ )0-|] + (y^+y2)cos 

Id-B,)"]! (31) 

★ TT 

In this section it is assumed that O<0^<^. Hence, if y^>vi 2 > the 
characteristic function is K-|(s) = 0 as given by (13), and if y^<y 2 
then in (13) y^ and y 2 are interchanged, where s-j is the relevant 
root with 0<s-|<l , 3' = l-s-j. 
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4. STRESS INTENSITY FACTORS 

For the crack tips fully imbedded in a homogeneous medium, 
(i.e., for a-|> 0 ,a 2 > 0 ) the stress intensity factors are defined by 

k3(b.) = lim^^b.l^2(r-b.)]TT.^^(r.6.) . 

k3(a^.) = , 62 = "^ + 9 ^ . 

i = 1 ,2 . (32) 

For the crack terminating at an interface, for example, for 
®1 ~ 0^ a 2>0 , k 3 is defined as 

l <3 = lim^^Q/? r^< (r,6^) . (33) 

0 Z 

For the crack crossing the interface, assuming that 0<6 i<tt/2 , 

at the intersection of the crack and the interface stress state 
is singular only on one side of the crack. In this case the stress 
intensity factor is ( arbi trari ly ) defined in terms of the inter- 
face shear stress as follows: 

k3(0) = lim^_^Q/2‘ r^f (r,Tr/2) , k3(0) = 0 , for y^>P 2 . 

0z 

k"(0) = lim^^g/F r ^1 (r,-V2) , k3(0) 0 , for yi<U2 

0 z 

(34a, b) 

In the numerical analysis the stress intensity factors are obtained 
in terms of asymptotic values of the functions f^(r) and f 2 (i^) O'" 
in terms of g.(+l) , (j = 1,2) (see equation (29)). 

J 

The calculated results given in the next section is based on 
constant crack surface tractions, namely 
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, -1 <x<l , i = 1,2. 


(35) 


q.(r)=Qi(x)= q.^ , a.<r<b. 

For example, if the nonhomogeneous plane is under uniform anti- 
plane shear 

' % < 36 ) 

then 

^10 " ’ ^20 " ^o^^"®l • 

Or, if the plane is under displacement loading 
9w-| 3 w2 

I3r = ay- " s 

then 

^IQ ~ ”^l^o^^^^l ’ *^20 ~ • (39) 


5. NUMERICAL RESULTS 

The numerical results for a crack lying in one half plane 
only (i.e., for a 2 = b 2 ) are given in Tables 3-5 where 

ag - (a-|+bi)/2 , c = a-j + Sg. (40) 

The tables show the normalized stress intensity factors defined 


by 


K( + l) 
K(-l) 




, for oag , 
for c = 3g 


(41) 


(42) 


The material pair of aluminum and epoxy is used as an example. 
In Table 3 Vi 2 ^Vi-] , in Table 4 F' 2 '^’^l ’ ^ Table 5 VI 2 “ ^ (the 
half plane). For the crack tip terminating at the interface 
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(i.e., for a-| = 0 or c = Hq) the corresponding power 6-| = 3 of 
the stress singularity is given in Table 1. 

The results for the problem of a crack crossing the boundary 
are given in Table 6. In this case the normalized stress intensity 
factors shown in the table are defined by 


KjCb,) 


l<3(b, ) 


Kj(b2> 


k,(b.) 


K-(0) 


k^O) 


(43a-c) 


In this example the external loads are assumed to be the displace- 
ment loading given by (38) and (39).. Note that since y-|<P2 
0<e-j<Tr/2 , k2(0) = 0. The related power of stress singularity 

8-| = 3' is given by Table 2. 

Figures 8 and 9 show some sample distributions of the density 
functions f^ and f2 for a-j = 82 = 0 which are normalized in the 
following form: 


Fi(t,) 




>1 3i+i 


-l<t^<l , 0<r<b^. , i = 1,2 (44) 

where is the load parameter defined by (38). Note that at r = 0 

or t. = -1 the density functions become unbounded having a common 
.3. 

r ' type singularity. 

Figures 10 and 11 give some idea about the crack opening dis- 
placement again for the case of a crack crossing the interface. 
Here, the normalized displacements W-j and W2 are defined by 

i ^1 3 i+7 . 

Wi-(t^.) = 2(-l)^ e^j(^) ^ [w|-w7J , i = 1,2 . (45) 
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23.08 for crack approaching 
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Table 6. Normalized stress intensities for crack through the interface, 
Uo/Wi = 23.08. 
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.0692 0.113 0.181 1.07 














It should be emphasized that ai the interface even though the 
displacements W] and W 2 are continuous, their derivatives are 
always discontinuous. This means that at r = 0 the crack 
opening displacement would always have a "kink". 
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Figure 3; Angular distribution of the measure of shear 

cleavage stress Tq in bonded half planes with 
a semi -i nfi n i te crack (see Figure 1 and equation 
5 for terminology ) . Material pair: Aluminum- 

Epoxy, 6 = it/8. 
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Figure 6. Angular distribution of the stress component 
for = tt/ 4 (see Figure 1 for terminology). 
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Figure 8. Di s tri butiorr of the normalized density functions 
in bonded half planes with a crack crossing the 
interface '(see equations 21, 2.9 and 44 and Figure 
7 for definitions). 




Figure 9. Same as Figure 8, ^ 2 ^^^ ~ 
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Figure 11. Same as Figure 10, ~ 


29 



